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Morphological transitions of phase separation associated with the asymmetry of lipid composition 
were investigated using micrometer-sized vesicles of lipid bilayers made from a lipid mixture. The 
complete macro-phase-separated morphology undergoes a transition to a micro-phase-separation-like 
morphology via a lorate morphology as a metastable state. The transition leads to the emergence 
of monodisperse nanosized domains through repeated domain scission events. Moreover, we have 
numerically confirmed the transitions using the time-dependent Ginzburg-Landau model describ¬ 
ing phase separation and the bending elastic membrane, which is quantitatively consistent with 
experimental results by fixing one free parameter. Our findings suggest that the local spontaneous 
curvature due to the asymmetric composition plays an essential role in the thermodynamic stabi¬ 
lization of micro-phase separation in lipid bilayers. 


Generally, binary liquids relax to complete phase sep¬ 
aration below a transition temperature. However, in a 
few systems, long-range repulsive interaction between mi¬ 
cro domains can lead to thermodynamic stable micro¬ 
phase separation Si- In terms of new material de¬ 
signs using membranes and the nanosized heterogeneities 
in cell plasma membranes M, one of the emerging 
problems is to elucidate how micro-phase separation in 
multi-component lipid bilayers can be stabilized and reg¬ 
ulated. One of the co-authors numerically predicts that 
the local coupling between the lipid composition and the 
membrane curvature can essentially lead to stable micro¬ 
phase separation [6[. Yanagisawa et al, have experimen¬ 
tally confirmed that the domain coarsening is transiently 
trapped due to the effective repulsive interaction between 
micro domains [3, [8|. However, to our knowledge, the 
mechanism of thermodynamic stable micro phase sepa¬ 
ration is not yet fully understood [1, S [sl • In this Let¬ 
ter, we demonstrate phase separation of a reconstituted 
two-component vesicle with an asymmetric lipid compo¬ 
sition between inner and outer layers, ubiquitous in cell 
membranes. Moreover, we numerically investigate the 
transitions using a dissipative dynamical model. 

Giant unilamellar vesicles (GUVs) were prepared ei¬ 
ther electroformation or natural swelling from a mixture 
of 35% DOPC, 35% DPPC and 30% Cholesterol at 55°C 
with pure water 0- Vesicles with excess area can re¬ 
sult in trapped coarsening after temperature quench 0. 
Thus, a small amount of Nal was added to the lipid film 
to generate vesicles with high sphericity ini- The average 
radius of the vesicles was 4.3 ± 2.0 jam (FIG. SI). Subse¬ 
quently, we decreased the temperature under the misci¬ 
bility transition temperature, leading to the spontaneous 
formation of liquid-ordered Lq domains (DPPC rich) 
in a liquid-disordered Ld matrix (DOPC rich), further 


* Email address; shimobayashi@chem.scphys.kyoto-u.ac.jp 
1 Email address; taniguch@cheme.kyoto-u.ac.jp 


coarsening to the fully phase-separated state through a 
collision-coalescence process. Simultaneously, GMl was 
dissolved in pure water and then ultrasonicated for two 
hours 0- A dynamic light scattering measurement indi¬ 
cates the formation of micelle structure (^ 6 nm) made 
of GMl molecules (^ 2 nm) in several tens of minutes 
after the ultrasonication (FIG. S2). Subsequently, 0.25 
mM DOPC/DPPC/Chol solution was mixed with an 
equal amount of 0.25 mM GMl solution. Time devel¬ 
opment of the average intensity of GMl was measured 
along the equatorial orbit of each GUV using confocal 
microscopy, indicating that a GMl insertion process into 
the GUV is equilibrated in approximately ten minutes 
(FIG. S3) [H, [ 23 , mi • Then, GMl preferentially parti¬ 
tions into Lo phase (FIG. S4). According to the intensity 
measurement of seven samples by the microscopy, the ra¬ 
tio of GMl concentration in Ld phase to that in Lq phase 
is estimated to be 0.49 ± 0.09. The average mole frac¬ 
tions of GMl in GUVs are estimated to be 2.90 ± 0.74% 
in Ld phase and 6.14 ± 1.70% in Lq phase, respectively, 
from comparisons between the mean intensities of GMl 
prepared in the film and inserted from the outside (FIG. 
S5). The following phase-separation transition was iden¬ 
tified by the difference of fluorescent intensities amongst 
the dyes used. The Ld phase, Lq phase, and GMl were 
stained by rhodamine-DOPE (0.5%), NBD-DPPE (1%), 
and Bodipy-GMl (0.6-1%), respectively 0- 

Typically, three different types of phase separation 
morphologies are observed: (a) a fully phase separated 
morphology (hereafter referring as “macro-phase”), (b) 
a lorate morphology (“lorate-phase”) and (c) a micro¬ 
phase-separation-like morphology (“micro-phase”). Lq 
domains in the lorate-phase and the micro-phase slightly 
bud toward the outside of the vesicles. At low Lq domain 
densities, they diffuse within a fluid Ld matrix (Movie 
SI and S2). On the other hand, at high densities, they 
organize into hexagonal order on the spherical surface, 
and hardly move (EIG. 1(a)). The time development of 
the fraction of vesicles exhibiting each one of these three 
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(a) macro-phase 





FIG. 1. Phase-separation transitions, (a) Confocal mi¬ 
croscope images of top (upper) and cross-sectional (lower) 
views of GUVs with three different morphologies (macro¬ 
phase, lorate-phase and micro-phase). Red and green rep¬ 
resent Ld and Lo phases, respectively. The scale bars are 10 
fim. (b) Time variation of the fractions of vesicles exhibiting 
each one of these three phases. The symbols denote the ex¬ 
perimental results and the solid lines are theoretically fitted 
ones. 


phases are shown in FIG. [Hb), where GUVs display¬ 
ing a coexistence of two different phases were discarded. 
All the vesicles initially exhibit a fully phase-separated 
state and subsequently the fraction of the macro phase 
state decreases continuously towards zero. On the other 
hand, the fractions of the lorate- and micro-phases in¬ 
crease from zero. The former reaches its maximum after 
about five hours and then decreases, but the latter con¬ 
tinues to increase. These results suggest that the macro¬ 
phase undergoes a transition to the micro-phase via the 
lorate-phase as a metastable state. Therefore, the succes¬ 
sive transition of state can be described by two coupled 
differential equations so that the experimental results 
are fitted as shown in FIG. 1(b) using the Levenberg- 
Marquardt method (Supplementary materials). The in¬ 
terface between and Lq phases are rapidly transformed 
at the transitions. The interface transformation at the 
transition from the macro- to lorate-phase could not be 
captured due to the fluorescence which suppresses it, pos¬ 
sibly due to the photo-induced oxidation of lipids |15j . 
On the other hand, we have succeeded in capturing the 


transformation from lorate- to micro-phase. A narrow 
constriction at a tip of a lorate domain emerges and im¬ 
mediately gives rise to domain scission (time scale less 
than 1 s; see white arrows in FIG. 2(a) and Movie S3). 
The scission event is repeated, which leads to the com¬ 
plete transformation from lorate domains to monodis- 
perse micro domains. The ratio of the length of the mi¬ 
nor axis of lorate domains, denoted by i?, to the radius of 
micro domains, r is the selected dimensionless wave num¬ 
ber. This was roughly estimated to be approximately 
0.7, which satisfy the Plateau’s prediction that the se¬ 
lected wave number is less than 1 (FIG.[2fb)). However, 
the instability should be related to not only the dissipa¬ 
tion of the interface energy, known as Rayleigh-Plateau 
instability, but also that of the bending energy of the 
membrane with non-zero curvature. Further investiga¬ 
tion of this instability represents intriguing future work. 
The size of the circular domains, 2r, were measured as a 
function of the Ld matrix radius, denoted by Ri (FIG. 
3). The solid line represents the best fit by a linear least- 
squares method, whose slope and intercept are 0.20 ± 
0.06 and 0.20 ± 0.17 /im, respectively. The results re¬ 
veal the emergence of nano-sized Lq domains near the 
limit of optical resolution and a linearly dependent rela¬ 
tionship between 2r and Ri , implying similarity of shape 
(FIG. 3). When a GMl solution of increased concentra¬ 
tion is added, we could not detect the micro domains, 
indicating that the domain size decreases below the op¬ 
tical limit. These findings strongly suggest that even a 
slight asymmetry of lipid distribution between inner and 
outer leaflets has the potential to dramatically change 
phase separation morphologies in lipid bilayers. 

Here, phase separation dynamics of an asymmetric 
two-component vesicle were numerically calculated us¬ 
ing purely dissipative dynamical model. The vesicle is 
represented by an one-dimensional closed line in a two- 
dimensional plane, parameterized by s. The order pa¬ 
rameter (j) is defined as the local difference between the 
concentration of DPPC, ^dppc, and that of DOPC, 
^DOPC, per unit area, i.e., <p = (/)dppc ~ ^dopc- Addi¬ 
tionally, for simplicity, we do not consider diffusion pro¬ 
cess of GMl inserted in the outer layer. In terms of 0, 
the two-component vesicle has the free energy functional 
U = -Fl + -F2: 


F,= J +PS, (1) 


F2 = 




( 2 ) 


Fi is the sum of the bending elastic and osmotic ener¬ 
gies where i<i{(j)) is the bending elastic modulus, H/2 the 
mean curvature and i7sp(0) the spontaneous curvature, 
p is the pressure difference across the vesicle and S the 
enclosed area. Note that the Gaussian curvature is zero 
because we assume translation symmetry parallel to the 
z-axis. We assume the following functional forms for /^(0) 
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FIG. 2. Domain scissions observed in the transition from 
lorate-phase to micro-phase, (a) The white arrows indicate 
narrow constrictions emerging immediately before each scis¬ 
sion event. The times are indicated for each image. The scale 
bar is 5 /xm. (b) Schematic illustration of domain scission. R 
and r represent the length of the minor axis of lorate domains 
and the radius of generated circular domains, respectively. 


and i?sp(0): + and i^sp(0) = 

where Hsp^ and Hsp^ are constants. F 2 is the 

Ginzburg-Landan free energy where /cp is the Boltzmann 
constant, T the temperature, a the unit cell length, <^o the 
correlation length. The potential /((/)) is a20^/2+a4(/)^/4, 
where a 2 and a4 are the negative and positive coefficients, 
respectively. 

We consider the time evolution of phase-separation 
transitions accompanied by dissipation of the free energy 
F. The position vector of a material point s at a time 
t is shown by r{s,t) = The simplest 

dissipative model reads 


dr{s, t) 
dt 


-r 


5r{s^ t) 


J Vgds'y{s,t) 


(3) 


where F is a kinetic coefficient and ^/gds the line 
element. The incompressibility condition, ^/g = 

y^{dx/ds)‘^ + {dy/dsY) = 1, for the local line element 
is guaranteed by a local Lagrange multiplier [6|. 

The following equations are derived from Eqs. o-©. 
See supplementary materials for the differential equa¬ 
tions for X and y. The equation for 7 ( 5 , t) is derived 
by substituting the differential equations of x and y in 
the relation = 0 as follows: 7 — — D = 0 with 

D = pH — xAy + yAx and Aa = Bd -b 2Bd where the 
dot on a variable X denotes dXjds and B is defined as 


(a) experiment simulation 
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FIG. 3. Size distribution of circular micro domains, (a) The 
definitions of the domain size, 2r, and the La matrix radius, 
i^i, in the experiment and simulation. The scale bar is 10 /xm. 
The domains where 0 > 0 (rich in DPPG) and (j) < 0 (rich in 
DOPG) are shown in green and red respectively, (b) Domain 
size, 2r, as a function of the La matrix radius, Ri. Both ex¬ 
periments (open circle) and simulations (solid circles) are in 
excellent agreement with the linear least-squares hts, which 
are the red solid-line and blue dashed-line respectively. The 
slopes and intercepts in experiments and simulations are 0.20 
±0.06 and 0.09 ± 0.18 /xm, respectively. In the simulations, 
the vesicle size is varied from 1.16 /xm to 4.56 /xm, corre¬ 
sponding to the change of N from 750 to 3000, for kq = 3000, 
Ki = 1/9, .ffip’ = 0.0765, and ffip’ = 0.0075. 


Secondly, we consider phase separation in lipid bilay¬ 
ers. The variation of ^(s, t) per unit time is given by the 
equation of continuity, i.e., d(j){sA)/dt = —{d/ds)J{sA) 
where J is the diffusion flux. The flux J is propor¬ 
tional to the lateral gradient of the local chemical poten¬ 
tial /x(s,t), namely J = —L(9/x(5, t)/9s), where L is the 
transport coefficient. The chemical potential /x is given 
by 5F/5(j){sA) and thus is calculated as follows: 


ld{s,t) = 


knT 


2 , dfi 1 Ok , 






(4) 

In our simulation, space, time, and energy are non- 
dimensionalized using = ^Q/|a2|, a^^/kBTL\a 2 \ 

and |a2|0eq^B^^^/<^5 respectively, where 0eq is de¬ 
fined as ^y\a 2 \/a 4 . This non-dimensionalization results 
in generating two scaled bending coefficients, kq = 
and Ki = /^i0eq//^O5 aud two scaled 

spontaneous curvatures, Hsp^ = Hsp^^ and Hsp^ = 
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(a) t = o t = ioo t = 500 



FIG. 4. Phase separation dynamics after setting a positive 
spontaneous curvature, (a) Cross sectional ID vesicles at 
100 and 500 in the phase-separation transitions. The domains 
where 0 > 0 (rich in DPPC) and 0 < 0 (rich in DOPC) are 
shown in green and red respectively, (b) Periodic prohle of 
0 as a function of the angle 6. The simulation parameters 
axe N = 2000, kq = 3000, ki = 1/9, Mp = 0.0765 and 
Mp = 0.0075. 

Hsp^cpeq^- We have attempted to use experimental re- 
suits for simulation parameters. The circular vesicle ra¬ 
dius at the initial state is estimated by A^^/27r, where N 
is a material point number. Thus, for example, by substi¬ 
tuting N = 2000 and a reasonable assumption ^ = 10 nm 
in this relation, the vesicle radius is obtained to be 3.2 
/rm, which is close to the average vesicle radius observed 
in the experiments (FIG. SI). The bending coefficient 
kq is difficult to determine from experimental values and 
it is the only undetermined parameter. On the other 
hand, the bending coefficient Hi is estimated to be 1/9 
by using the ratio of the bending coefficient in DPPC 
rich phase to that in DOPC rich phase, 1.25 [T6|. The 
spontaneous curvature of the bilayer in each phase is ob¬ 
tained by the relationship where 

and are the spontaneous curvatures of the in- 
ner and outer monolayers p/7|, 118| . We assume that the 
spontaneous curvature of a monolayer is the linear func¬ 
tion of the mole fractions of components naii^. Using 
the reported spontaneous curvatures of DPPC, DOPC 
and GMl (0.003 A“^ dH, -0.006 [H and 0.017 
A“ 0, respectively) and the measured mole fractions 
of GMl in two phases (^ 6% in Lq phase and ^ 3% 
in Ld phase), the spontaneous curvatures of the bilay¬ 
ers in two phases are obtained. Using the above values, 
Hsp^ and Hsp^ are estimated to be 0.0765 and 0.0075. In 
our model, we have ignored the flip-flop phenomena of 
GMl having a bulky and hydrophilic head group as even 
smaller lipids without such a head group show very slow 
reaction rates (several hours to a day) . 

In our simulation, we set fully phase-separated circu¬ 
lar vesicle with zero spontaneous curvature, i.e., Hsp^ = 0 


and i7sp^ = 0. The vesicle was sufficiently relaxed to an 
equilibrium state. At t = 0 we assume that GMl was 
instantly inserted into the outer layer, which results in a 
change of the spontaneous curvature in each domain. We 
set the non-zero spontaneous curvatures, Hsp^ = 0.0765 
and = 0.0075, and observed the following dissipa¬ 
tive process from t = 0 to t = 500. Note that the instant 
insertion of GMl is reasonable as the time scale of in¬ 
sertion is negligible compared to that of the subsequent 
phase separation processes. Based on the framework, we 
performed numerical simulations using the following pa¬ 
rameters: N = 2000, Ho = 3000, Hi = 1/9, Hsp^ = 0.0765 
and = 0.0075; so that we observed the gradual phase 
separation transitions with shape deformations (FIG. 4, 
Movie S4). In FIG. 4(a), the domains where 0 > 0 (rich 
in DPPC) and (p < 0 (rich in DOPC) are shown in green 
and red respectively. The monodisperse Lq domains are 
distributed at approximately regular intervals along the 
line and slightly bud toward the outside of the vesicle, 
which are in excellent agreement with experimental re¬ 
sults (FIG. 4). Qualitatively speaking, the transition 
results in the increase of the boundary energy F 2 , but 
also the decrease of the bending elastic energy, Fi, as 
the budding of Lq phase causes the membrane curva¬ 
ture to become a closer value to the spontaneous curva¬ 
ture, therefore total free energy decreased. Also, note 
that the bending energy is decreased more by budding of 
the stiffer Lq domains both in experiment and numerical 
simulation. Moreover, we compared Ri and 2r obtained 
from our experimental results with those from our numer¬ 
ical calculations. The value Ri has been defined as the 
average length between the center of the circular vesicle 
and each material point which satisfies p < 0 (FIG. 3(a)). 
The value 2r has been defined as the average length be¬ 
tween two neighboring material points at which the sign 
of (p is changed as a function of 0 (FIG. 3(a)). The an¬ 
alyzed values and those standard deviations are plotted 
in FIG.3(b) with the best fitting (the blue dashed-line), 
where the material point number alone has been varied 
from N = 750 to N = 3000 for Hq = 3000, Hi = 1/9, 
Hsp^ = 0.0765, and = 0.0075. Quantitatively, the 
numerical results are mostly consistent with experimen¬ 
tal results. Furthermore, we investigated a domain num¬ 
ber A'cap, ‘2r/Ri and total free energy as a function of 
the mole fraction of GMl in Lq phase. The ratio of the 
mole fraction in Lq phase to that in Ld for N = 1000, 
Ho = 3000 and Hi = 1/9, which is kept constant, equals to 
two. We have found that there are discontinuous jumps 
of Acap and 2t jRi as a function of the mole fraction. Al¬ 
though Acap = 1 below the transition fraction, indicat¬ 
ing a vesicle of macro-phase, Acap increases and 2r/Ri 
decreases as the mole fraction of GMl increases above 
it (FIG. S6). However, the detail of the discontinuous 
transition remains poorly understood. 

To summarize, we have experimentally found that fully 
phase-separated phase with asymmetric lipid composi¬ 
tion undergoes a transition to the micro phase via the 
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lorate-phase as a metastable state, resulting in the emer¬ 
gence of monodisperse submicron scale micro domains. 
Numerical analysis based on the bending elastic model 
and TDGL equation shed some light on the mechanism 
of the transitions, implying the significance of the spon¬ 
taneous curvature due to the asymmetric composition 
in the thermodynamic stabilization of micro-phase sep¬ 


aration in lipid bilayers. Our findings could give some 
important clues as to the physical principle for new ma¬ 
terial designs and behind the heterogeneities existing in 
cell membranes. 
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